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* Abstract. We have investigated the role of plasma instabilities in the origin of a solar flare. 

Three types of stability problems are discussed: the magnetohydrodynamic instability of an 
infinitely conducting fluid, the instability due to the magnetic field gradients, and the finite 
conductivity instability. These instabilities are examined in a situation that occurs when two 
plasmas trapped in a pair of sunspot magnetic fields approach each other. The assumption of 
infinite conductivity gives rise to stable situations. The instability growth time in the magnetic 
field gradients that exist near sunspots is extremely large. The growth time for finite conduc- 
tivity instability is of the order of a few seconds or minutes, and we therefore believe that the 
flare can be produced by finite conductivity instabilities in the solar atmosphere above a 
group of sunspots. Our analysis does not help explain the acceleration of charged particles to 
high energies; for this, the problem of acceleration must be considered separateiv. _ 
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1. Introduction 

It has been suggested [ Parker , 1962; Dungey, 
1961; Gold and Hoyle , 1960] that instabilities 
may play an important role in explaining the 
origin of solar flares. Dungey [1961] has sug- 
gested that a pinch instability can take place in 
the ionized solar atmosphere near the neutral 
point of the sunspot magnetic fields and also 
that runaway electrons may develop in the 
transition region, leading to a two-stream in- 
stability. The possibility of the pinch instability 
is based on the assumption that a current may 
be driven along the neutral point (or line) be- 
tween two regions of oppositely directed mag- 
netic fields to such an intensity that the self- 
magnetic field of the current overpowers all 
other forces. Although this may be a possibility, 
our analysis in section 2 shows that such an in- 
stability does not arise in the equilibriums con- 
sidered here. Gold and Hoyle [I960] consider a 
particular geometry of a complex sunspot group 
(in relation to the problem of origin of flares) 
in which a pair of bipolar sunspots is parallel to 
another pair of bipolar sunspots. We shall con- 
sider here the simpler geometry of four sunspots 
each of which belongs to a pair of conjugate 
spots arranged in a straight line. This geometry 
is considered following Sweet [1958]. 

Severny [1958, 1959] and Severny and Shab- 
hanski [1961] have worked out a model of the 
pinch effect in the neighborhood of the neutral 
point of the sunspot fields assuming that strong 


currents are generated by the magnetic field 
gradients. Severny arrives at the conclusion that 
the plasma being compressed by the two sun- 
spot fields is unstable when the magnetic pres- 
sure is greater than the kinetic pressure of the 
plasma. The result is that if a compression be- 
gins it will proceed at an ever-increasing rate; 

1. e., the pinch effect takes place. The phenomena 
that take place in the chromospheric or coronal 
plasma near the sunspot region are quite dif- 
ferent from the pinch effect produced in the 
laboratory. Figure 1 shows an equilibrium con- 
figuration of the pinched plasma produced in 
the laboratory. Figure 2 shows one of the pos- 
sible equilibrium configurations that can arise 
in the solar plasma before the occurrence of a 
flare; this case is considered in detail in section 

2. The basic difference between the two cases is 
that, whereas in the first the magnetic field out- 
side the plasma is a decreasing function of the 
distance, in the second it is uniform. We shall 
show in section 2 that the equilibrium configura- 
tion considered in Figure 2 is stable. The phys- 
ical reason for the instability of the configura- 
tion in Figure 1 is that any plasma moving away 
from the axis of the cylinder is at a point at 
which the magnetic field, and therefore the mag- 
netic force, is less, and the perturbation con- 
tinues to grow. Similarly, if any part of the 
plasma moves toward the axis of the cylinder it 
is further pushed by the magnetic field due to 
the increased magnetic force outside the plasma 
cylinder. If we now consider the equilibrium 
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Fig. 1. Pinched plasma produced in the labora- 
tory. 

configuration of Figure 2 we at once see that 
any displacement of the plasma in or out of its 
region is in the presence of the same magnetic 
force and there is no cause of instability as in 
laboratory pinches. The same holds for the 
equilibrium configuration sketched in Figure 3. 

In sections 3 and 4 we shall analyze some of 
the known instabilities that apply in the present 
situation, with the understanding that an in- 
stability with growth rates of the right order 
constitutes the flare phenomena, as an insta- 



Fig. 2. A possible equilibrium configuration be- 
fore the onset of flare. The magnetic field on both 
sides of a zero field region is oppositely directed. 
Outside is assumed vacuum. The plasma, the cur- 
rent sheath, and the magnetic field extend both 
ways, perpendicular to the plane of the paper, up 
to infinity. The current in the sheath flows in the 
Y direction. 


bility helps to blow the material up. The dis- 
cussion presented here does not help in under- 
standing how the particles get accelerated to 
energies of the order of a few million electron 
volts; for this, acceleration processes must be 
considered separately. 

2. Magneto hydrodynamic Stability of 
Equilibrium Configurations 

Observations show [Saverny, 1958] that, before 
the occurrence of a flare, the sunspots are 
approaching each other with the velocity of 
a few kilometers per second. Assuming that we 
have a group of four sunspots arranged in a 
straight line as considered by Sweet [1958, 
Figure 3], let us suppose that a magnetohydro- 
dynamic equilibrium is reached when the two 
magnetic fields of the sunspot groups are sepa- 
rated by a distance 2x 0 with a region of zero 
magnetic field and hot plasma between. For 
simplicity, we neglect the gradient of the mag- 
netic field on both sides of the zero field region. 
The possibility of such an equilibrium con- 
figuration will be increased if it is proved that 
the magnetic field gradients do not give rise to 
any serious instabilities. In section 3 we will 
show that the e folding growth time of any 
perturbation due to the magnetic field gradient 
instabilities is of the order 10 10 ’ 3 seconds. This 
allows us to neglect any instabilities due to the 
gradients of the magnetic field. The pressure 
anisotropy instabilities do not arise in this case 



Fig. 3. Another possible hydromagnetic equilib- 
rium of oppositely directed magnetic field re- 
gions. The plasma region and the current sheath 
extended to infinity in Y, Z directions. Right-hand 
half-plane is filled with plasma with magnetic field 
(0, 0, Bo) in cartesian coordinates; the left half is 
filled with plasma with oppositely directed mag- 
netic field (0, 0, — Bi). 
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because the collision frequency is sufficiently 
high in the chromosphere as well as in the 
corona to preserve the pressure isotropy during 
the compression of the plasma. For coronal 
plasma of N ~ 10 8 particles/ cm*, T ~ 10 6 °K, 
the frequency of electron-electron collisions is 
of the order of 2.5 X 10* per second and the 
mean free path of an electron is of the order of 
1 km. Since the time involved in the flare process 
is of the order of minutes, and the length is of 
the order of thousands of kilometers, the pressure 
will be assumed isotropic throughout the dis- 
cussion. 

The basic equations in this section are the 
magnetohydrodynamic equations as used by 
Kruskai and Schwarzschild [1954] and as further 
employed by Jaggi [1962], With the isotropic 
pressure the equations of motion are 


grad p = 

- j X B 

c 

(1) 

<*P _ 
dt 

“ di v (pv) 

(2) 

dB 

dt 

curl (v X B) 

(3) 

1 dp 
p dt 

ydp 
p dt 

(4) 

di v B = 

0 

(5) 

curl B = 

(47r/c)j 

(6) 


where the perturbations of the physical entities 
in the plasma are of the form (e.g., the perturba- 
tion in the pressure denoted by subscript 1) 

Pi = Pl (x)e ik ' v+ik ''- i<at (10) 

Linearizing equation 4 gives 

pi = f pi ( ii ) 

where f = V ypo/po is the velocity of sound in 
the plasma, and p Q and p 0 are the equilibrium 
pressure and density of the plasma, respectively. 
We easily find that 


B x = 0 (12) 

pi satisfies the equation 

^5 — (fci 2 + &2 2 — «*/£*) Pi = o (13) 
and has for its solution 

Pi = Ae Kx + Be~ Kx (14) 

with K = V + k 2 2 (o3 2 /? 2 ). 

Also 

Pi = f(Ae Kx + Be~ Kx ) 
v z = ( Ae Rx - Be~ Kx ) 

upo 

v, = — {Ae Kx + Bf Kz ) > 


(15) 


(16) 


where the displacement current has been neg- 
lected. As a possible equilibrium configuration 
for the plasma let us assume that there is no 
magnetic field in the plasma confined between 
two parallel planes x = For x > z 0 there is 
vacuum and the magnetic field (0, 0, B 0 ) is 
constant and parallel to the surface. Taking 
the z axis parallel to this magnetic field, we have 
a current in the y direction at x = +x 0 . For 
x < x 0 the magnetic field is (0, 0, — B 0 ). The 
boundary conditions are 

n- [B] = 0 (7) 

Ip + (£*/ Sir)] = 0 (8) 

If we perturb the boundary of the plasma- 
vacuum interface as 


x = + i 


(v^)x »i a w + »'*»«— 

* e 

co 


‘ (9) 


v. = — (Ae Kx + Be~ Kt ) 

w Po 


The solutions for perturbation in magnetic 
field in the vacuum x > x Q and x < x 0 are, 
respectively, 


and 


Bi. 


DKye' K ' x 

B\y 

= 

Dik 2 e~ K ' x 

Bi, 

= 

Dik 3 e~ K ' x 

Bi, 

= 

D l K 1 e KlX 

B\y 

= 

D,ik 2 e K,x 

B im 

= 

D 1 ik 3 e K,x . 


(17-18) 


where K x = \ / 'k 2 2 + and D and D x are 
constants of integration; the factor »+«,«-<»! 
has been omitted on both sides of equations 14 
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and 17-18 for convenience. Applying the boundary 
conditions (7) and (8) at x = +x 0 and x — — x 0 , 
we obtain the dispersion equations 


or 

2 

CO 


K x 


V A % 2 k 

K x 

V/h 2 k 


(e 2Kx ° 


- 1) = (e Kx ° + 1) 
+ 1) = (e 2Kr ° - 1) 


(19) 


which show that oo 2 cannot be negative, and 
therefore the equilibrium is stable. 

The other equilibrium configuration considered 
by us is shown in Figure 3. For x > 0, we have 
a uniform plasma of density p 0 , pressure p 0 , and 
the magnetic field (0, 0, Bo), For a: < 0, we 
have a uniform plasma of density pi, pressure pi, 
and magnetic field (0, 0, — Bi). At x — 0, we 
have a current (c/47r)(£i — Bo) in the y direction. 
Following the analysis given above, we obtain 
the dispersion equation in the form 


k 2 2 + h 2 + 


k 2 + k 2 + 



co 2 co 2 



1/2 


where 

Co = VtPo/Po 

V 0 = B 0 / V4 irpo 



Ci = VtPi/Pi 
7, = BJ V47TP, 


If co 2 < 0, the left-hand side is positive and 
the right-hand side negative. This shows that 
no solution is possible for oj 2 < 0; i.e., no in- 
stabilities exist. For B 0 = B h p 0 = pi, p Q — p i, 
i.e., the symmetrical case, the dispersion equation 
is of the form 


k 2 k 2 




2 

co 

7 2 2 
& 3 C 0 



- 0 ( 21 ) 


which gives a quadratic equation having two 
positive roots. Thus there are no instabilities of 
the current layer in equilibrium between two 
ionized clouds having oppositely directed uniform 
magnetic fields. If we take the gradient of the 
magnetic field into account we predict that the 
current layer will still be stable because the 
magnetic field on both sides has increasing 
magnitude. 


3. Magnetic Field Gradient Instabilities 

We shall now show that there are no serious 
magnetic field gradient instabilities for long 
wavelengths to affect the assumption of the 
equilibrium configurations studied in the previous 
section. Krall and Rosenbluth [1962] have con- 
sidered a plasma in equilibrium in a magnetic 
field of the form 

B = £«( 1 + «)i, (22) 

where i t is a unit vector in the z direction. This 
magnetic field configuration is quite suitable for 
our problem. The distribution function for the 
electrons and ions can be assumed to be 


/ \Ir.V /2 

= n( -))( — ) exp [—a,(vf + »„*)]• 

exp [—^i{v,)]g (23) 


where the subscript j refers to electrons ( j = e) 
and ions (j — i ) ; g satisfies the equation 


dg _ eB 0 
dx me 


[1 + ex] = 0 

dvij 


(24) 


Choosing 

g.= l- «‘(* + hx* ~ f) 

!7, = 1 - «'(» + W + f) 


where e 1 = ef? 0 2 /[4xtt&(7\- + T e )], n being the 
number density of the particles, S2 e — eBo/mc , 
0^ = eB Q /Mc , m = the mass of the electron and 
M that of an ion, Krall and Rosenbluth arrive 
at unstable solutions of the form e Bt with the 
growth rate s given by 


R. («) 



{[EiiXTi/T,)} 2 + x 2 } 2 


*- 




( 25 ) 
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for k R. < 1, t 1 / ( 1, and 

, 0 m /2 ( 2xT > \ f (i 3 xt '\ 

4n, Vr,/ P V T . ) lFl \2’2’ T. ) 


ft. w = 




(26) 


for k R 0 > 1, €]/ € 1, where R e and R { denote 

the Larmor radii of electrons and ions respec- 
tively. Here X = 2d i $l i 2 /k e should not be 
confused with the coordinate x. k is the wave 
number of purely longitudinal electrostatic oscil- 
lations of the form 


E - (a/| k | )e iU ' x ’E i (XT i /T e ) 

~ e XTJT ’/{XT,/T.) 

for large values of X, and iFi is a confluent 
hypergeometric function. For the case of sunspot 
magnetic field gradients, 


e 


3 dB 
B dx 


1 

10 1( 


B 0 ~ 10 3 gauss 
X ~ 0 .243(FV7\) 


where X = 2 t r/k. For T = 5 X 10 3 °K, X ~ 1 
cm, X ^ 10 13 , i.e., an extremely large quantity. 
Equation 25 yields 


R. (s) 


Hh 21 10 36 


10 " 


(27) 


i.e., an extremely small number. 

Conversely, if k R e > 1, i.e., wavelengths less 
than the electron cyclotron radius, the growth 
rate is extremely large, of the order of inverse of 
the growth rate for k R, < 1. This shows that 
only extremely small wavelengths are unstable. 
The Debye length in the lower chromosphere is 
(kT/4ime 2 ) lli ~ 5 X 10~* cm, and the electron 
cyclotron radius corresponding to 5 X 10 3 °K 
temperature is of the order of 10“ 3 cm. Any 
wavelengths less than the electron cyclotron 
radius will also be less than the Debye length, 
and the theory may not apply. (This is satisfied 
in the corona as well.) At any rate, the instability 
of extremely small wavelength, although taking 
place instantaneously, will not be efficient to 
enhance the thermal and electrical diffusion of 
the plasma. 


The magnetic field gradient instabilities dis- 
cussed so far in this section are only important 
much before the equilibrium between the mag- 
netic pressure and the kinetic pressure is reached 
(if it is reached at all) because of the condition 
€ l /e ^>> 1 imposed on solutions 25 and 26. After 
the compression continues, the temperature of 
the plasma builds up; the pressure P begins to 
increase until it is comparable with the magnetic 
pressure £ 2 /8t. The kinetic pressure develops 
more rapidly than the magnetic field pressure 
[Severny, 1958, p. 321], and an equilibrium 
configuration can be reached when P is com- 
parable with B 2 / 87r. 

It might be thought that the pressure ani- 
sotropy instabilities can develop in this case, 
but it can be shown that owing to the slow rate 
of approach of the two sunspots and owing to 
the high collision frequency the pressure ani- 
sotropy cannot develop to the extent of pro- 
ducing an instability. As has been mentioned 
earlier, the collision frequency in the lower 
corona is of the order of about 10* per second, 
and, since the compression time is of the order 
of hundreds of seconds, the collisions will keep 
the pressure quite isotropic in the present 
situation. It is only if the phenomena take place 
in about 10*“* second that the pressure anisotropy 
can be important. To calculate the magnetic 
field, the density, and the temperature after the 
compression we may use the adiabatic equations 
deduced by Chew , Goldberger , and Low [1956]. 
With P fl = jP ± , these equations are 

PB 2 , . 

= constant (28) 

p 

— constant (29) 

1\ £> 

from which we obtain 

£ 3 

-^2 = constant (30) 

If L 0 , Bo t Po denote the initial length, magnetic 
field density, and the pressure of the plasma, 
the quantities B x N x P x for a given L x can be 
calculated from NiL x = N 0 L Ql B x 3 /N x 2 — £ 0 */X 0 f , 
Pi/Ni 61 * = Po/Nf* For B 0 = 10 3 gauss, 

No = 10 l2 /cm 3 , L 0 = 10 10 cm, T 0 = 5 X 10 3 °K, 
we obtain for L x = 10 7 cm, N x = 10 18 /cm 3 , 
Bi = 10 s gauss, Ti = 5 X 10 5 °K. The collision 
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frequency remains about the same as before the 
compression. 

4. Resistive Instabilities 
The instabilities discussed so far apply in the 
case of infinite conductivity of the plasma under 
consideration. Recently some work has been 
done assuming finite conductivity of the plasma 
[Furth et al . , 1963]. We shall use the results to 
make an estimate of the growth rates of resistive 
instabilities. 

The infinitely thin current layer as sketched 
in Figure 3 will now be assumed to have finite 
thickness denoted by a. We will also assume a 
magnetic field in the current layer. With finite 
conductivity, equation 3 has to be replaced by 

^ = curl (v X B) — curl curl bJ (31) 

where y is the isotropic resistivity of the plasma. 
The magnetic field in the current sheath has 
the form 

B - i.B a (z) (32) 

Further assuming the current layer to be incom- 
pressible (this amounts to using div v = 0 in 
place of equation 2 and deleting equation 4), 
two unstable modes are deduced by Furth, 
Killeen, and Rosenbluth for waves propagating 
in the z direction. They are called the rippling 
mode and the tearing mode. The rippling mode 
arises in the case in which the resistivity varies 
with distance in the current layer; the growth 
rate of the instability which is the inverse of the 
e folding time of a perturbation is given by 

R. (s) = a 2/5 £ 2/5 / 4ra\r (33) 

where a (— l/rj) is the conductivity of the 
plasma in emu 

m2 

S = i0“ 2 y 1 

V 47T P o/ Bq 


T in electron volts, a = ka , where k is the wave 
number of waves in the current layer (propa- 
gating in the current layer). Also this instability 
occurs for S~ 2n < a < S 2/3 . Using a = 2 X 10“ 14 
T*'* emu, T = 10* °K, a = 10* cm (10* °K ^ 86 
ev), we have S = 7 X 10 3 , and the e folding 
time which is the inverse of R € (s) is of the 
order of 10 3 /& 2/3 . Now, from S~ in < a < S 21 * 
we get 2.8 < l/k < 10*. For long wavelengths, 
the e folding time for the rippling mode is large. 
Only for small wavelengths, of the order of a 
few centimeters, the times are reasonable and 
the instability of rippling mode is important for 
the problem in hand. 

The tearing mode is unstable for a < 1 or 
wavelengths greater than the thickness of the 
current sheath. In this case the e folding time is 
of the order 10* k 2/6 , where k > 1/2500. This 
time will be small for long wavelengths, and the 
order of magnitude is also suitable for the 
problem. 

The above results are derived on the assump- 
tion that the current layer is incompressible, 
and it is proved by Furth, Killeen, and Rosen- 
bluth that the results are not much different if 
compressibility is taken into account. 

To arrive at the above figures, we have used 
the temperature of 10* °K in the higher chromo- 
spheric region. Some people [ De Jager, 1959, 
1961] are of the opinion that the temperature 
of the chromosphere region does not rise very 
much and that only the density becomes higher. 
A differing view, based on theoretical models, is 
adopted in section 3 of this article, following 
Severny [1958]. For this reason we have evalu- 
ated the growth rates of rippling and tearing 
modes in Table 1. The magnetic field at the 
neutral point, where the current sheath is 
formed, will be weak. This is another reason 
why the growth rate can drop off by several 
orders of magnitude. 


TABLE 1. Rippling and Tearing Mode Instability 
If t r denotes the e folding time for the rippling mode, this table evaluates t r (k) 2,& = 2.96 X 10 -2 T 9I10 /B Q 2/6 . 
The e folding time for the tearing mode is obtained from t t =t T 10 4 (fc) 4/5 . 


T, °K 

B = 10 
gauss 

B = 10 2 
gauss 

B = 10 3 
gauss 

B = 10 4 
gauss 

B = 10 5 
gauss 

B = 10* 
gauss 

10 4 

44.4 

18.6 

7.42 

2.96 

1.17 

0.469 

10 6 

3.7 X 10 2 

1.48 X 10 2 

59 

23.5 

9.35 

3.72 

10* 

2.96 X 10 3 

1.17 X 10 3 

4.69 X 10 2 

1.78 X 10 2 

74.03 

29.6 

10 7 

3.72 X 10 3 

1.47 X 10 3 

5.9 X 10 2 

2.24 X 10 2 

93.02 

36.58 
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For the rippling mode the e folding time is 
obtained by multiplying each member of the table 

by (k-'y*.- We find that ™ > (ir') 2,s > 

for the whole range of the table. The upper 
numbers apply for B — 10, T = 10 4 ; the lower, 
for B = 10 6 , T = 10 7 . An intermediate value of 
1 shows that suitable time scales are obtained 
for almost all the values of B and T in the table. 
The above results are deduced for N ~ 10 15 
particles/ cm 3 . If, instead, we use N ~ 10 12 
particles/cm 3 , it amounts to multiplying each 
number in the table by 10 3/8 ~ 4 . N occurs in 
t r or t t as A r ~ 1/5 . 

For the tearing mode the e folding time is 
obtained by multiplying each member of the 
table by lO 4 ® 2 ' 8 , where, for the whole range of 

the table, 10~ 2 > (&) 2/5 > The upper 

number applies for B = 10, T ~ 10 4 ; the lower, 
for B = 10®, T = 10 7 . Assuming k 2l& to be of 
the order of 10 -3 , we see that the e folding times, 
for large magnetic fields and small temperatures, 
are reasonable. We therefore conclude that the 
resistive instabilities give reasonable times for 
the problem in hand. 
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